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using a soft cascade approach. We derive an evolution equation for the 
gluon density in soft gluons cascades in the leading logarithmic appro- 
ximation of perturbative QCD, and show that this equation reproduces 
BFKL equation in the forward case. The whole cascade is equivalent 
to a single gluon whose self-energy is responsible for gluon reggeization. 
The same type of equation is obtained for the QED case. 
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1 Introduction 

Investigation of hard scattering amplitude in the kinematics where invariant 
energy s is large while transfered momentum t is fixed is an important problem 
of QCD. The solution of this problem in the Leading Logarithmic Approxima- 
tion (LLA), collecting terms of the form {as In s)" , was proposed about twenty 
five years ago, in the BFKL equation jTj. These original papers were based on 
an effective triple gluons vertex and a bootstrap condition in t— channel. This 
approach shew the relation between perturbative QCD and reggeon calculus 
which was proposed a decade before. The main feature of this work is the 
reggeization of gluon, which appears not to be elementary but composite, as 
being a pole in complex momentum plane, with color octet quantum num- 
ber. Another important result is that Pomeron occurs as a bound state of two 
reggeized gluons in singlet color channel. Although BFKL equation looks like 
a ladder type equation, it effectively sums up an infinite number of t— channel 
gluons. 

A few years ago a dipole picture was proposed [21 El IH El El E! • It is based 
on subsequent emissions of soft gluons, which iterates in the large Nc limit. In 
this approach elementary degrees of freedom are made of colored dipole. The 
equation for dipole density in this soft cascade turns out to be identical to the 
BFKL equation, although it was found in the multicolor limit, while original 
BFKL was derived for finite Nc- 

The aim of this paper is to show that BFKL dynamics can be simply re- 
covered in the forward limit only using this classical soft emission vertex used 
in the dipole model jSj. We deal with gluons at finite Nc rather than with 
dipoles at large A''^. We obtain BFKL equation as an evolution equation for the 
gluon density in the cascade. This equation is similar to DGLAP pH ITI H ITT ] IT^ 
except that the evolution variable is Inl/x rather than InQ^ . Some related 
approaches were elaborated in jT3| and [T^ . 

Our paper goes as follows. After preliminary recalling linear classical cas- 
cades properties, we generalize them in order to incorporate non abelian effects 
in the LLA approximation. The main idea is that despite the complicated in- 
ternal structure of the cascade, the amplitude of soft emission is the same as 
the one off a single gluon. Using this idea we rederive BFKL equation for 
forward case. The same technique, when applied to QED, gives the same type 
of equation. 

2 Classical cascade 

Consider the scattering of two particles. The asymptotic behavior of the cross 
section is determined by the parton cascade from incoming particles or, in other 
words, by their wavefunction on the light cone. For large invariant energy s 
(the so-called Regge limit), the amplitude is dominated by t-channel gluons 
rather than quarks, because of their largest spin. As a consequence gluons 



dominate in every channel, and quark contributions will be neglected below. 

We shall consider a semi-hard region where gluons can be treated to be soft 
with respect to the large invariant energy s, although still perturbative. It is 
known that in any massless field theory scattering amplitudes are dominated 
by IR contributions of two types, namely collinear and soft singularities. This 
semi-hard region reveals soft gluons contributions, which sum up into LLA. 

The physics behind soft resummation is that the field of an ultrarelativistic 
source can be treated as a state of almost real particles while the vertex for 
real or virtual soft gluon emission can be taken to be the same. 

Let us denote pa_ and pb the momenta of colliding particles A and B. 
Neglecting their masses in the high energy limit we can take p\ = p'b = so 
that the invariant energy reads s = 2pA ■ Pb- Vectors pa, Pb can be used for 
the Sudakov decomposition of any vector as k = apA + Pps + k±. 

Each particle develops a parton cascade, so the scattering of the incoming 
particles A and B could be in principle reduced to elementary parton processes. 
However the complete description of the cascade is a very complicated pro- 
blem involving all perturbative as well as nonperturbative effects. There are 
certain limits, where the cascade looks more simple. A well known example 
is deep inelastic scattering. For large virtuality Q^ the parton density obeys 
the DGLAP evolution equation collecting the powers of as InQ^. Although 
predictions based on the DGLAP evolution yield good agreement with the 
present experimental data, this approach fails for parametrically small Bjorken 
variables x, when the powers of 0:5 In 1/x dominate. This is the region of 
Regge kinematics where the total invariant energy is much greater than the 
other invariants, including the virtuality of the deeply virtual photon. The 
leading In 1/x behavior of the hard amplitude is given in this region by the 
BFKL theory. 

An important feature of the Regge kinematics is that the partons are soft 
there in the sense that the main contribution appears for small longitudinal 
variables, a,/3 -C 1. Only soft emissions from the external incoming or out- 
going particles should be taken into account in this approximation. Since 
the momenta of all incident external particles are supposed to be along pA or 
Pb directions, two peculiar gauges are of special interest, namely the gauges 
Pa ■ A = and pb ■ A = 0. These gauges suppress soft emission from pa and 
Pb lines respectively. In the following we shall deal with emission from pb line 
in the gauge pa ■ A = 0. 

The vertex for emitting a soft particle is rather universal and determined 
by the classical current proportional to the momentum of the source. It can 
be easily obtained, for example, from the triple gluon vertex in the soft limit. 
Consider emission off the incoming particle B, as illustrated in Fig.H 

This incoming line is attached to an amplitude (shown as a blob in Fig.P) 
whose peculiar form plays no role so long as we are dealing with soft emission. 
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Figure 1: Soft vertex for incoming line. The blob symbolizes any amplitude 
connected to the incoming line. 



The soft vertex reads 



r^(fc) 



PBje\k) 
PbK — 10 



Pb ■ £\k) ^ 
9 — s ^^c 



(1) 



where e^{k) is the helicity vector of the soft emitted or absorbed gluon carrying 
momentum k and color index c. The matrix T'^ depends on the representation 
of the color group. For gluon like object it is expressed through the structure 
constants, T^,„ = ifa'ca, the indices c, a', a being in the adjoint representation. 



Cn; "'n 



C2, ^2 






Ci, ki 



Pb, a 



Figure 2: Amplitude for n gluon emission or absorption off an incoming gluon- 
like object, represented as two thin lines. 



The amplitude for emission or absorption of n gluons (see Fig.|21) is given 
by the product of the elementary vertices. 
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Here the field A is the asymptotic free field 
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whose positive and negative frequency part corresponds respectively to emis- 
sion and absorption of gluons. Using the relation 
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we get 

= — j— j dzi--- dZnP[nB ■ A{znnB) ■■■ ub ■ A{zinB)]^,a , 
with 

and where the symbol P means the conventional ordering of the fields. The 
sum over the arbitrary number of soft gluons is given by a P-exponent along 
the momentum direction of the parent particle: 
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A similar computation for outgoing line results in replacing in the previous 
formula j_^ with /q°° . These formulae are nothing more than the Wilson line 
taken on the trajectory along the momentum of the incoming (resp. outgoing) 
particle. 

3 Inclusion of interaction among emitted par- 
ticles 

The previous computation only takes into account emissions from the parent 
line which is the incoming line (resp. outgoing). The above formulae are strictly 
speaking only valid for QED. They do not take into account the interaction 
between emitted particles. In QCD, in the soft approximation, these correc- 
tions correspond to subsequent decays of the emitted gluons. Indeed in this 
approximation only end lines emissions contribute, and thus to get the total 
cascade tree-level amplitude each field in the P-exponent has to be replaced 
by the P-exponent itself. Generally this turns into a complicated non-linear 
equation. Another problem is to incorporate loop corrections. This second 
problem will be considered in section 0] 

The first problem, the partons subsequent decays, simplifies in the Regge 
kinematics implying the longitudinal momenta to be small in the formula Q, 
a,/3 <C 1, while 

as, Ps ^ k']_, ajSs ~ fc^ , (4) 

where /c^ is a typical transverse momentum scale. This kinematics ensures the 
soft vertex to be of the form (^ regardless from where a particular gluon is 
emitted off or where it is absorbed in. Indeed, the emission of a gluon with 



momentuin /c2 = (^Pa + (3pb + ^2± off the parent gluon carrying momentum 
ki = apA + bpB + ki± is again given by a soft vertex similar to (^ 

A., N ^ fci • e\k2) ^ 
^cik2)=9^^:^T,. (5) 

In the hght cone gauge pa ■ A{x) = 0, the polarization vector reads 

^iik) = e\{k) ■ PA, + e\,{k) (6) 

with 

The emitted gluons are on-shell. Anyway, the gluons which dominate the 
amplitude in the Regge limit are soft and thus quasi-real, therefore the on- 
shell condition is unessential. The transversality condition for these quasi-real 
gluons reads fc ■ e = 0, which implies for the helicity vector (jH)) 

e\{k) = 2^-^^. (7) 



In this soft gluon approximation, (3 <^ b and assuming the typical transverse 



momenta to be of the same order, fc? i ~ A;, i , the numerator of the soft vertex 



is approximated as 



ki ■ e\k2) =he\-- ki^ei ^ bps ■ e\k2). 

Using the mass shell conditions, a = kfj_/bs, a = k"^^/ l3s, the denominator 
also simplifies as 

s s 

ki-k2 = [ajS + ba)- - kij_ ■ k2± ^ ba- . 

It follows that the vertex (jH)) can be written as 

T^M) = /-^^^T^, (8) 
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which has exactly the universal form (0). It confirms the physical picture that 
soft emission is determined by the total current of the source rather than by 
its internal structure. The vector pb is the momentum of the source. 

With this universal vertex combined with Jacobi identity for the color ma- 
trices T*^ the emission of the new soft particle from the "ends" a and c in the 
Fig. 1 looks like if it were effectively emitted off the line a' at the left from the 
particle c. This is illustrated in Fig. El 
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Figure 3: Illustration of the soft gluon universality using Jacobi identity. 

Since the emitted gluon is real, /c^ = a/3s — fc^ = 0, and thus, combining 
Eqs.® and (|7j), the vertex (jHJ can be rewritten in the form 



T^M = 9T 



2ks_ ■ e± 



(9) 



Repeating the new emissions from all " ends" in the same manner we come back 
to the Fig. 121 where now the gluons are ordered according to their (3 value, the 
smaller (3 being on the left of the largest ones. 
Note that the decreasing order of /3 variables. 



A > /?2 > ■ ■ ■ > /?„, 
implies in the Regge kinematics increasing order of a's, 

tti ^ ^2 ^ ■ • • ^ «n, 
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so that ai + a2 — 0^2, ai + - ■ ■ + an — an- This property allows to present the re- 
sulting soft cascade tree amplitude through an ordered exponent in momentum 
space, 

^tree = Pp exp \J dp d'^k± V{p, k±)\, 



with 
and 
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In the gauge pa- A = 0, this equation reduces to 
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where we have used expression Q- This expression incorporates both emission 
and absorption of the partons described by the hght cone creation, a^^(/5, q±), 
and annihilation, axc{P,q±) = a^^{P,q±) = a|'^(— /?, — gj_), operators labelled 
by the longitudinal momentum fraction (3 > 0, transverse momentum q±, po- 
larization and color indices A and c, which satisfies 

[aAc(/5,g±),a+,,(/5',gl)] = 2P 6xy6,,, 6i(3 - /3')6^'\q^ - q'^). 

The ordering symbol P/j means that the fields A{i3) with the smallest /3 
value are on the left from the fields with the largest ones. The minimal value 
P = x\ plays the role of an infrared cut-off in the amplitude ^tree- Note that 
in the case of incoming line, the variable "time" z in the P-exponent (jH)) is 
Fourier conjugated to the longitudinal Sudakov variable a, and thus ordered 
oppositely. Since on the other hand, for mass-shell particles a ~ 1//5, it turns 
out that \z\ ~ (3. Thus, P-exponent defined with respect to z is in accordance 
with Pp. 

When the cut-off value is lowered, that is when we take x'^ < x\, it allows 
for the emission or absorption of an extra soft particle whose longitudinal 
momentum lies in the interval x\ < (3 < x\. This yields the new amplitude 
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the amphtude ^tree standing as the source for the new soft particles. The 
whole tree amphtude can be symbolically presented as an infinite product of 
elementary emissions or absorptions in the infinitesimal intervals Ax, 



^tree = IT 1 + T d(3 fd^k^V{P, k, 

^ L Jx—Ax J 



(15) 



4 Virtual corrections and evolution of the cas- 
cade wave function 

The previous form is convenient to include the virtual contributions. Besides 
the amplitude to emit or absorb a real gluon at each elementary step there 
is an amplitude for the case when the new phase volume in the longitudinal 
space remains non-occupied. In other words it means that the emitted gluon is 
reabsorbed at the same step, which corresponds to a loop correction. It does not 
changes the number of particles, but it is responsible for the renormalization 
of the whole state. 

Despite the composite internal structure of the whole system of cascade 
and source, within the soft emission approach this system looks like the source 
carrying given momentum and color. The virtual correction results into self- 
energy insertion into the source propagator. In particular, if the source is gluon 
like, then the full system looks like a gluon. In this case, the virtual correction 



results into self-energy insertion in the gluon propagator. In LLA and for the 
cascade plus source system of momentum g ~ Pb + 5'±, the propagator of this 
system is modified as 
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with the function tt{xx, qY) determining z-factor, or normalization of the state. 
To keep the value of the norm fixed, the cascade plus source wavefunction is 
multiplied by z,^^{q±), where q±_ is the total transverse momentum of the 
system whose constituents fill the interval of longitudinal momenta [xa, 1], 
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The contribution of the soft gluons appearing at the next evolution step in 

-| Icy 

the interval [x\ — 6xx, xx] should be compensated by the factor z. '_g^ ^ , {q±)- 
Basically for a finite interval one can write a loop expansion of the form 
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only the first variation determining the evolution equation, 
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In LLA function uj{q±) is given by a one loop diagram calculated in the axial 
gauge. 

Before discussing an explicit form of the virtual contribution note that the 
momentum q± in the argument of uj is the total transverse momentum of the 
system of cascade and source, which emit a soft particle as a whole. Introducing 
the total momentum operator Pi, uj{P±) obviously gives u!{q±) when acting on 
a whole system of cascade plus source state of total momentum q±. Thus, in a 
similar way as in the tree case (see Eq. (jl4p ). one step of evolution reads 
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The full cascade S matrix can be symbolically written through the elemen- 
tary steps product^ 



1 r 
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X Jx-Ax 



(17) 



^Here Sxx denotes the variation of the cut-ofF value while Aa; is taken as a notation for an 
infinitisemal step in the infinite product representation of the cascade S matrix. Principally 
one can put Sxx — Aa;. 



where the brakets are ordered from the smallest /3's values at the left to the 
greatest at the right. The operator 

'^^(^x) = n[i-^(i^±)— + r dp fd^k^v+ip, k^) 

^, L X Jx—Ax J 

acting on the vacuum creates the soft constituents of the cascade, while the 
operator 

picks out the components with a fixed transverse momentum. Let us in- 
troduce the operator b^^{pB, q±) for creating the bare source, having trans- 
verse momentum q_j_, helicity and color indices a, a, and denote this state as 
\pB,q±,a,a)s = 6+„(pb, g±)|0). Then the state 

\PB,Xx,q±,a,a) = d^p±^+{xx, pi_)b^^^{pB, q± - p±)\0) 



describes the system of cascade plus source with given total transverse mo- 
mentum q±, cut-off value xx, and source helicity and color. The total tranverse 
momentum of the cascade is given by p±. The operator P± acting on the right 
results at every step into the transverse momentum of the total system at the 
previous cut-off value. 

Changing xx we evidently have for the variation of the cascade plus source 
state 

0X\ 

S\pB,xx,q±,a,a) = LL!{q±)\pB, q±, xx, a,a) (18) 

2, 



+ / dp d k±V^{P,k±)\pB, (q- k)±, Xx, a,a). 

xx-Sxx 

Let Q be an operator, which probes the parton distribution in the cascade, 
for instance,^ 

Q = — '^^/±/(a;,/±)a+(x,/_L)a^c(a;,/±) (19) 

•'XQ Zx 

with a weight function f{x,l±). The average value is generally expressed 
through the density function nj, depending on the cascade total transverse 
momentum and cut-off, 

{pB,Xx,q±,a,a\Q\pB,Xx,q'±,a',a') = Sa,a'5a,cT'5^'^\qi. - q'_i)nf{xx,q±). 



^In principle the operator Q could probe both the cascade content and the source, without 
modifying the following discussion. One could add to Q a term 

1/2 J d^l^f{h)b+JpB,l±)b,,c{PB,l±) 

acting on the source. If the source is gluon-like this is the same as to include in the weight 
function /(cc, l±) in Ea. H19|l a term proportional to 5(1 — x). 



Its evolution with x\ value is determined by the variation of the state ()18|). If 
we suppose that the operator Q has a natural cut-off xq ^ xa, then the new 
emitted gluon does not affect the operator vertex (the operators ax^dP, k±) 
commute with Q for xx — 5xx < (3 < xx) and the variation of the matrix 
element decays into the sum 

Auf = AiUf + A2n/, 

as illustrated in Fig.EJ The first term is given by the extra particle amplitude 





c^^+^ ^c>+o 




Figure 4: Variation of the density with respect to the cut-off, in the case 
Xq 3> x\. The double thick line stands for the system of source plus gluon 
cascade. 

squared times the average of the operator Q over the rest part of the amplitude 
corresponding to the previous cut-off value xx and recoil transverse momentum. 
It corresponds to the first diagram in Fig.|3]and can be written symbolically 

Ain/(xA,g±) = {pB,xx,q±- k±,a,a\ax,cQat^c\PB,xx,q'±- k±,a,a') (20) 

or in explicit form 

Ainf{xx,q±) = 2Nc-——^ ^^^ J -p-'^/(^A' ^^ ~ ^^)- (21) 

Note that the previous equation is written in the case where the source is in 
the adjoint representation. In the general case, Nc should be replaced by the 
Casimir of the corresponding representation. 

It is important to note that in soft loops calculations gluons can be consid- 
ered as real, in accordance with the fact that there is almost no difference be- 
tween real and virtual massless soft particles. Indeed, consider the loop where 
one single gluon is emitted and reabsorbed, as illustrated in the Fig. HI The a 
integral which occur in the loops can be closed around the pole of the emitted 
gluon propagator. For, using the same technique which leads to the eikonal type 
expression ((T)), the denominator of the integrand reads (a — i(5)^(«/3s — k\ + i5) 
for the first graph of Fig. 0] or (a — i5){al3s — k'j_ + i6) for the second and third 
graphs. Both denominators leave the singularities in a plane on the opposite 
side of the real axis. The numerator of the gluon propagator in light-cone 
gauge reads 

4A;2 
d^^u{k) = -el{k)el{k) - -^pAf^PAu- 
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When performing the a integral, if among the two poles one choose to close 
around the physical pole of the gluon, the second term drops out since it cancels 
the a singularity in the propagator. Thus, only the first term remains. One 
then immediately gets the soft universal vertex squared Q for the first graph. 
The second term in the variation arises due to virtual correction, 

A2n/(xA,g±) = -21n4cu(gx). (22) 

It is illustrated by the second and third diagrams of Fig. HI 

Differentiating with respect x\ we arrive at the evolution equation 

xx-^nf{xx,q±) = "^A/'cT^^y -^—^^nf{xx,k±) - 2uj{q±)nf{xx,q±) ■ 

In the case where xq is smaller than the typical value of xx, there appears 
an additional term when the extra soft gluon operator is contracted with the 
operator Q, as shown by the additional fourth diagram in Fig.El It gives the 





C^^+ ^ ^ o+o 





Figure 5: Variation of the density with respect to the cut-off, in the case 
Xq < Xx- 

following contribution to the variation of Uf : 

2 

Asnf{xx,q^)=2N,-f—fixx,q±)\n^. (23) 

[Z TC) X^ 

The full inhomogeneous equation thus reads 



2 j2 1 

-2N,-^J -^—^Ufixx^k^) - 2uiq^)njixx,q±). (24) 

If one is interested in the density number of gluon of momentum r± , the weight 
function / in Eq. (fTU|) should be chosen as 

f{x,h) = l6'{h-r^), (25) 

the 1/2 factor being related to the average over gluon transverse polarization. 
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^^dx. 




+ 




Figure 6: BFKL equation for the gluon density in a cascade. The double hne 
stands for the source. 



In that case, Eq, 
d 



turns into the BFKL equation (for t = 0) 



xx 



dx) 



G{xx, r±, q±) = - Nc 



(2 7r)3 ri 



-2Nr 



9^ 



£k^ 



(27r)3y {q-k)\ 



■G{xx, r_L, k±) - 2uj{q±)G{xx, r±, q±) . (26) 



The inhomogeneous term, which is of lowest order, can be interpreted as the 
initial gluon contribution in the absence of cascade. Since we consider here 
gluons emitted from a source, it contains a factor N^g^ /{2tiY, multipliying 
1/r^ which is reminiscent of the two t— channel gluon propagators (one being 
compensated by polarization contribution after angular averaging of fc^). 

The equation Eq. (PUJ) looks like a Bethe-Salpeter equation with a kernel K 
acting from below, which is illustrated in Fig.El 

Actually the density G is only a function of p^ — r± due to translational 
invariance. It means that Eq.(|^ could equally be written as a Bethe-Salpeter 
equation with a kernel acting from above, by just an exchange between of 
variables r±_ and p±. 

A natural physical value for the IR cut-off xx is xx ~ fi'^/s where /i^ is some 
typical scale for tranverse momentum. This follows from the Regge kinematics 

©• 

Using Eqs. lfT^ and (jHSJ, one can finally obtain the unintegrated parton 

density from G as 



g{xx,r±, q± 
d 



- — {pb, xx,q±,a,a\ - X1"<^,c«^,c|Pb, Xx,q'±,a',a' ) 



X,c 



dx) 



G{xx, r^, q_i) 



(27) 



In this formula, r± is the tranverse momentum at which parton density is 
measured while q± is the total tranverse momentum of the cascade plus source. 
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5 Computation of uj 

5.1 Evaluation of uj based on the gluon cascade univer- 
sality 

The expression (OT|l provides a simple way to find the function u}{q±) based 
on universahty. From the point of view of emitted gluon, the whole system 
corresponding to the previous cut-off value plays the role of the source, and 
the 2;-factor is prescribed to this whole system. On the other hand, this whole 
system looks like a gluon (in the case where bare source is gluon like). Thus, 
the z-factor for the whole system and for a single gluon should have the same 
functional form. What we need is the virtual correction for the whole system, 
but we shall calculate instead the z-factor of the dressed emitted gluon. The 
inclusion of ^-factor corresponds to the replacement 0^,0- — > z^^"^ ciRx,a in the 
expression (PU]) . With z-factor included the formula (PT|l reads 

A,„,,fe.,J = 2iV.j|^l„-l^ (28) 

/dp 
-^-Z[x^-5x^,x^]{pL)nf{xx,q±-p±). 
P± 

On the other hand there is a dual description of the same dressed gluon as a 
composite cascade state spread over the interval [x\ — 6x\, x\] and carrying 
transverse momentum pj_. The two particles component of this state is given 
by the amplitude to emit two soft gluons off the source (which is the whole 
system corresponding to the previous cut-off value), 

dx f dp f d'^k± k± ■ 6^2 



2(3 J (27r)3/2 kl 

xx-Sxx xx~5xx 



I X \ . r ■ r [ dx r dij r 

IPb, g±, xx -dxx, aa)2 = ^gJac^d^gJdc2b J ^ J ^ J 

X (-27r)3/2 ^ ^ (- - t)2 "' ^^'"^ ^^' ^^^ ^^^^^ ^^' P±-k±)\pB, q±-p±, xx, b,a), 
the gluons being ordered with respect to their longitudinal momenta, 

xx — Sxx < P < X < Xx, 

which is reflected in the integration limits. 

Averaging the operator with this amplitude we arrive at the following ex- 
pression for the gluon density variation 

A,n,{xx,g^) = 2A^^(|^^l-'^;^/^-/(^A,gx-px) 



(27r)3 J " klip - kf_ 
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Note that the double integration in x and /3 results into 1/2 In^ — ^ — . This 



xx-Sxx ' 



factor 1/2 reflects the Bose symmetry of the two gluons system. Recalling 
Ea.dTHl) and Eq.(I2HI) we get 



2uj{p_i 



Nr 



(2vr)= 



d'ki 



vl 



k±{p — k) 



2 • 



(29) 



5.2 Direct calculation 

The formula (j^^ can be compared with a direct calculation of the polarization 
operator. The exact one loop result in the dimensional regularization {D = 
2 + 2e) and axial gauge is 

.2 / D^ rl 



7r(p±) 
where 



-2N, 



9 

c o 2 



r 1 



YJio^^l 



l-(3)pi/An]--\l-(3)K{f3), (30) 



K{P) 



P 



+ ^ + (3{l-(3) 



1-/3 p 

is the DGLAP kernel. Gluon self-energy is known to coincide in the axial gauge 
with a Sudakov form factor whose double logarithmic behavior originates from 
the product of transverse and longitudinal divergencies. The latter one comes 
about when one gluon momentum in the loop becomes soft, /3 — »• 0. To separate 
the longitudinal logarithms needed in order to find the function uj{p^) we keep 
only the singular piece in K{P) and cut the integral in (jHOI) at /? = x\, using 
dimensional regularization only for tranverse divergency, which turns into 

-^ ' D\ , ^ . D, .1 , . . 1 



T^iixx, P_l) 



-2N, 



9 



c o 2 



r 1 



J ) r±/4vr 



— -1 



In — = 2lu(p±) In — . 

Xx Xx 



(31) 



This results into 

2uj{p^) 



2Nr. 



87r2 



1 , pi 

- + In ^ + 7ij 



while dimensionally regularized integral ()29|) gives 



2u{qi 



Nr. 



11 



[ql/ATi 



— -In 



2-£ 

2 



2Nr 



87r2 



- + In^ +7s 
e 47r 



+ 0(e) 



r^(f-i) 
r(L) - 2) 

+ 0(6). 



Thus, although functions uj{pi_) defined by eq. (j^^ and eq. (j^ look different 
their nonvanishing parts are equal. 

Further, if we consider the renormalization of a heavy color source of mass 
M due to classical current emission, the divergent part of the source self energy 
amounts to the same expression 



t^m{xx,p\) 



-2N, 



9' 



57r" 



r 1 



with source mass M and virtuality M^ 
or composite nature of the source. 



n 
2 

n2 



M^ 



P 



An 



In 



Xx 



(32) 



p ~ pj_ regardless of an elementary 
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6 QED cascade 

The above treatment can be applied to electrodynamics. Generally speaking, 
the wave function of the source, which could be any charged object, for instance 
an electron, is made of the bare source itself surrounded by a cloud of soft 
photons. The operator Q can probe either the source or photon component of 
the wave function. If we introduce a cut-off for the longitudinal momentum of 
the soft photons, we can study the effect of changing this cut-off as we did for 
the QCD case. We get the same three type of contributions as illustrated in 
Fig.l^l where the wavy lines now represent photons. 

The z-f actor is related in this case to the renormalization of the source, 
and can be obtained from Eq. ()H2|l . One needs only to replace g by the electric 
charge Z e and put iVc = 1 . The real contribution is only due to the emission 
off the source line, since photon has no charge. In the corresponding Wilson 
line Q, the ordering of the fields is not important since they commute at equal 
light cone time. It is related to the fact that in QED, the soft cascade structure 
does not assume any additional (3 ordering (fTTH) . 

The equation satisfied by the corresponding nj is the same as Eq. ()24j] . after 
performing the replacement g'^Nc — > Z^ e^ . In the case of photon density, one 
gets 

Xx-^—G^{xx, r±, q±} = -77^-^ 2 

~^727r)3 J ( -kY '^^^^^^ ^^' ^^^ ~ '^^^^^'> '^7(3^A, r±, q±) . (33) 

This equation has just the same form than the BFKL equation (f^H|l for QCD. 
Let us stress once more that contrarily to the QCD case, uj{q±) corresponds 
to the reggeization of the source and not of the photon, which is of different 
nature. In QCD, in the case where the source is gluon like, uj{q^) can be 
associated either to the source or to the cascade. We have used this property 
to get the expression for uj{q±). From the point of view of source dressing, it is 
obtained through the one loop gluon self energy (see Eq.ijHH)). Equivalently, it 
is obtained through the normalisation of the wave function (see Eq.(|29|)). 



7 Cascade wave function 

Let us now turn back to the amplitude (J17p . Despite the symbolic character 
of the infinite ordered product it can be presented in a closed form. Con- 
sider to this end the amplitude to emit n quasi-real gluons with momenta 
^n, Q±i, ■ ■ ■ ,Xi, gj_^, Xi > X2 > . . . > Xn- It is given by the matrix element 

■'- A„,c 



n 5 ^n 5 



,Xl,cv,a',cT;a',a{Xn,q±n, ■ ■ ■ ,Xl,q±i; g_|_, ^J.) (34) 

( I K',a'{PB, q'±) aA„,c„(a;„, g±J ■ ■ ■ ax,,c^{xu g±J ^^{xx) K^^{pb, g±) I ) . 
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Expanding the brackets, the product in the $"'"(xa) operator can be reorganized 
as 



«^^(:^a) = n 

X 

X n 



1 - UJ{P) 



Ax' 



X 



XV— Aa; 



d(3n d^k^„V+{(3n,k^J 



^n— 1 



1 - ^(P) 



Ax' 



X 



^2 r -^ At 

n i-^(^) 



Xl 



X 



r d/3i fd^k^,V+i/3,,k^JY[ 

Jxi—Ax J J. 



^ ■" - Ax' 



X 



+ 



The terms which are not exphcitly written gives zero contribution to the matrix 
element flH4|l when performing the contractions with a operators. Using the 
fact that 



""' ■" -^ Ax 



n 



X 



exp \ — I -^uj[F) 



X2 P 



and that the operator P resuhs into the total momentum of the state occuring 
to the right, we get for the amplitude 



An; Cn, ..., A 



l,ci;a',a\a',a\^l: Q'Xl; • • • ) ^ny Q±n] Q^i Q±) 



a;^\ ^(g±„+---+'?±i+9x) 






^(qi 



^g±i+<?x) 



Xrj, / 



Ul, 



, X-n—l 



re::,H?±„-J 



xr^;(?±.; 



^iq±) 



{ I b^'^a'iPB, q'±) btJpB, g±) I ). 



This formula, in the LLA approximation, provides all information on the cas- 
cade wave function. It resembles multi-reggeon formula but it includes only 
soft vertex. 

8 Conclusion 

This paper is based on the ordering of longitudinal variables in soft cascades, 
which means energy ordering of emitted particles. It differs from coUinear 
approximation which implies angular ordering jT^] . In this sense, this treatment 
describes the dynamics responsible for BFKL evolution, while the collinear 
one is related to the conventional partonic DGLAP picture. We have shown 
explicitly that the infrared evolution equation for the parton density in the soft 
cascade reproduces forward BFKL equation. 

In both BFKL and DGLAP equations the evolution is written with respect 
to a cut-off. DGLAP involves a transverse cut-off Q^ due to collinear singular- 
ities, while in our case the cut-off is longitudinal, due to soft singularities. 

In DGLAP case the unintegrated parton distributions, which depend both 
on the longitudinal and tranverse momenta, are related to the derivative of the 
structure functions with respect to the transverse cut-off. Similarly we have to 
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consider the derivative with respect to x\ of the function G(xx,r±,p±), which 
satisfies BFKL equation, in order to get the unintegrated distribution. 

Our approach is based on a kind of duahty in the description of a gluon 
at high energy. From the viewpoint of the soft vertex the gluon state can be 
treated as a single particle with given momentum and color. On the other 
hand the gluon has internal structure, a cloud of many soft gluons surrounding 
it. It looks like a composite object with a nontrivial wavefunction. Hence a 
rather simple interpretation of reggeization as a manifestation of soft cascades 
structure appears in this picture. The trajectory ui{q±) is nothing more than 
the gluon Z factor in the axial gauge. Both in QCD and QED, this function is 
associated with the charged source self energy. In DGLAP case the divergencies 
occurring at x — > 1 are regulated by the l/(x — 1)+ prescription, corresponding 
to UV Z-factors due to UV renormalization of parton wave functions. In BFKL 
case the IR divergencies for k± —>■ disappear because of virtual corrections 
which can be treated as IR renormalization of the gluon wave function. The 
equation (j^ can be wieved as describing 2 reggeized gluons, or, equivalently, 
2 gluon soft cascades. 

The extension of this approach for non forward case and for Generalized 
Leading Log Approximation (fHJEI^l) will be presented in a forthcoming 
paper. It would be also interesting to study more complicated structures like 
for example the triple Pomeron vertex. 
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